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Abstract 

In many natural synchronization phenomena, communication between individual elements 
occurs not directly, but rather through the environment. One of these instances is bacterial 
quorum sensing, where bacteria release signaling molecules in the environment which in turn 
are sensed and used for population coordination. Extending this motivation to a general non- 
linear dynamical system context, this paper analyzes synchronization phenomena in networks 
where communication and coupling between nodes are mediated by shared dynamical quan- 
tities, typically provided by the nodes' environment. Our model includes the case when the 
dynamics of the shared variables themselves cannot be neglected or indeed play a central part. 
Applications to examples from systems biology illustrate the approach. 

Keywords: Synchronization, quorum-sensing, systems biology 

1 Introduction 



Many dynamical phenomena in biology involve some form of synchronization. Synchronization 
has attracted much research both from the theoretical, see e.g. (Strogatz 2003), (You et al. 2004), 
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(McMillen et al. 2002) to cite just a few, and experimental (Yagamuchi et al. 2003), (Pye 1969) 
viewpoints. The particular case of synchronized time-periodic processes, where time-scales can 
range from a few milliseconds to several years (Winfree 2001, Newman et al. 2006), includes e.g. 
circadian rhythms in mammals (Gonze et al. 2005), the cell cycle (Tyson et al. 2002), spiking 
neurons (Izhikevich 2006) and respiratory oscillations (Henson 2004). 

When modelling such networks, it is often assumed that each node communicates directly with 
other nodes in the network, see e.g. (Park et al. 2008, Bohn & Gracia-Ojalvo 2008) and references 
therein. In many natural instances, however, network nodes do not communicate directly, but 
rather by means of noisy and continuously changing environments. Bacteria, for instance, produce, 
release and sense signaling molecules (so-called autoinducers) which can diffuse in the environment 
and are used for population coordination. This mechanism, known as quorum sensing (Miller & 
Bassler 2001, Nardelli et al. 2008, Ng & Bassler 2009) is believed to play a key role in bacterial 
infection, as well as e.g. in bioluminescence and biofilm formation (Anetzberger et al. 2009), (Nadell 
et al. 2008). In a neuronal context, a mechanism similar to that of quorum sensing may involve 
local field potentials, which may play an important role in the synchronization of clusters of neurons, 
(Pesaran et al. 2002, Boustani et al. 2009, Tabareau et al. 2010, Anastassiou et al. 2010). 

From a network dynamics viewpoint, the key characteristic of quorum sensing-like mechanisms 
lies in the fact that communication between nodes (e.g. bacteria) occurs by means of a shared 
quantity (e.g. autoinducer concentration). Furthermore, the production and degradation rates of 
such a quantity are affected by all the nodes of the network. Therefore, a detailed model of such 
a mechanism needs to keep track of the temporal evolution of the shared quantity, resulting in an 
additional set of ordinary differential equations. 

Mathematical work on such quorum sensing topologies is relatively sparse (e.g., (Garcia-Ojalvo 
et al. 2004, Tabareau et al. 2010, Russo & di Bernardo 2009b, Katriel 2008)) compared to that on 
diffusive topologies, and often neglects quorum variable dynamics or the dynamics of the environ- 
ment. This sparsity of results is somewhat surprising given that, besides its biological pervasiveness, 
quorum sensing may also be viewed as an astute "computational" tool. Specifically, use of a shared 
variable in effect significantly reduces the number of links required to achieve a given connectiv- 
ity (Tabareau et al. 2010). 

In this paper, we derive sufficient conditions for the coordination of nodes communicating 
through dynamical quorum sensing mechanisms. These results can be used both to analyse natural 
networks, and to guide design of communication mechanisms in synthetic or partially synthetic 
networks. We first consider, in Section IXT| the case where the network nodes (e.g., the biological 
entities populating the environment) are all identical or nearly identical. We then focus, in Section 
13. 2[ on networks composed of heterogeneous nodes, i.e., nodes of possibly diverse dynamics. In this 
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case we provide sufficient conditions ensuring that all the network nodes sharing the same dynam- 
ics converge to a common behavior, a particular instance of so-called concurrent synchronization 
(Pham & Slotine 2007). In Section |33| the results are further extended to a distributed version of 
quorum sensing, where multiple groups of possibly heterogeneous nodes communicate by means of 
multiple media. Finally, in Section |U we propose a strategy for controlling the common asymptotic 
evolution of the network nodes. Section [S] studies the dependence of synchronization properties 
on the number of nodes, a question of interest e.g. in the context of cell proliferation. Section E] 
illustrates the general approach with a set of examples. 

Our proofs are based on nonlinear contraction theory ((Lohmiller 8z Slotine 1998)), a viewpoint 
on incremental stability which we briefly review in Section |51 and which has emerged as a powerful 
tool in applications ranging from Lagrangian mechanics to network control. Historically, ideas 
closely related to contraction can be traced back to (Hartman 1961) and even to (Lewis 1949) (see 
also (Pavlov et al. 2004, Angeli 2002), and e.g. (Lohmiller & Slotine 2005) for a more exhaustive list 
of related references). As pointed out in (Lohmiller & Slotine 1998), contraction is preserved through 
a large variety of systems combinations, and in particular it represents a natural tool for the study 
and design of nonlinear state observers, and by extension, of synchronization mechanisms (Wang & 
Slotine 2005). 

2 Contraction theory tools 
2.1 Basic results 

The basic result of nonlinear contraction analysis (Lohmiller & Slotine 1998) which we shall use in 
this paper can be stated as follows. 

Theorem 1 (Contraction). Consider the m- dimensional deterministic system 



where f is a smooth nonlinear function. The system is said to be contracting if any two trajectories, 
starting from different initial conditions, converge exponentially to each other. A sufficient condition 
for a system to be contracting is that there exists a constant invertible matrix B such that the so- 
called generalized Jacobian 



X = f{x,t) 



(1) 




(2) 



verifies 



3A>0, Va;, Vt > 0, fi{F{x,t)) < -X 
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where /i is one the the standard matrix measures in Table 1. The scalar A defines the contraction 
rate of the system. 

For convenience, in this paper we will also say that a function f{x, t) is contracting if the system 
X = f{x, t) satisfies the sufficient condition above. Similarly, we will then say that the corresponding 
Jacobian matrix ^{x,t) is contracting. 



Table 1: Standard Matrix measures 



vector norm, ■ 




induced matrix measure, ^ {A) 




1 1 v~^m 1 
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We shall also use the following two properties of contracting systems, whose proofs can be found 
in (Lohmiller & Slotine 1998, Slotine 2003). 

Hierarchies of contracting systems Assume that the Jacobian of ([1]) is in the form 

(3) 

corresponding to a hierarchical dynamic structure. The Ja may be of different dimensions. Then, a 
sufficient condition for the system to be contracting is that (i) the Jacobians Jn, J22 are contracting 
(possibly with different G's and for different matrix measures), and (ii) the matrix J12 is bounded. 

Periodic inputs Consider the system 

x = f{x,r{t)) (4) 

where the input vector r{t) is periodic, of period T. Assume that the system is contracting (i.e., 
that the Jacobian matrix |^(x,r(t)) is contracting for any r{t)). Then the system state x{t) tends 
exponentially towards a periodic state of period T. 

2.2 Partial Contraction 



dx 



x,t) 



Jll J12 
J22 



A simple yet powerful extension to nonlinear contraction theory is the concept of partial contrac- 
tion (Wang & Slotine 2005). 
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Theorem 2 (Partial contraction). Consider a smooth nonlinear m- dimensional system of the form 
X = f{x,x,t) and assume that the so-called system y = f{y,x,t) is contracting with respect to 
y. If a particular solution of the auxiliary y-system verifies a smooth specific property, then all 
trajectories of the original x-system verify this property exponentially. The original system is said 
to be partially contracting. 

Indeed, the virtual y-system has two particular solutions, namely y{t) = x [t) for all t > and 
the particular solution with the specific property. Since all trajectories of the y-system converge 
exponentially to a single trajectory, this implies that x (t) verifies the specific property exponentially. 

2.3 Networks of contracting nodes 

This section introduces preliminary results on concurrent synchronization of networks, which will 
be used in the rest of the paper. 

Consider a network consisting of heterogeneous nodes: 



where Ni denotes the set of neighbors of node i and 7 is a function defined between two set of 
indices (not necessarily a permutation), i.e. 



Thus, two nodes of ([5]), e.g. Xi and Xj, share the same dynamics and belong to the p-th group 
(denoted with Qp), G Qp, if and only if 7 (i) =7 (j) = p. The dimension of the nodes' 

state variables belonging to group p is n^(i), i.e. Xi G R"t(') for any Xi G Qp. In what follows we 
assume that the Jacobian of the coupling functions are diagonal matrices with nonnegative 
diagonal elements. We will derive conditions ensuring concurrent synchronization of i.e. all 
nodes belonging to the same group exhibit the same regime behavior. 

In what follows the following standard assumption (see (Pham & Slotine 2007) and refer- 
ences therein) is made on the interconnections between the agents belonging to different groups, 
(Golubitsky et al. 2005). 

Definition 1. Let i and j be two nodes of a group Gp, and if they receive their input from elements 
i' , j' respectively, then: (ii) i! and j' belong to the same group Gp>; (ii) the coupling functions 
between i-i' and j-j' are the same; (Hi) the inputs to i and j coming from different groups are the 
same. If these assumptions are satisfied, then nodes i and j are said to be input-equivalent. 




(5) 



7:{l,...,iV}^{l,...,s} s<N 



(6) 
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Given this definition, we can state the following theorem, which generalizes results in (Pham & 
Slotine 2007) to the case of arbitrary norms. Its proof is provided in the Appendix. 

Theorem 3. Assume that in ^ the nodes belonging to the same group are all input- equivalent and 
that the nodes dynamics are all contracting. Then, all node trajectories sharing the same dynamics 
converge towards each other, i.e. for any Xi, Xj G Qp, p = 1, . . . ,s, 

\xj (t) — — )■ as t — +oo 

In the case of networks of identical nodes dynamics, the above result amounts to only requiring 
contraction for each node. 



3 Main Results 

In this Section we present our main results. We first provide sufficient conditions for the syn- 
chronization of a network composed by nodes communicating over a common medium, which is 
characterized by some nonlinear dynamics. We then extend the analysis to a number of cases, by 
providing sufficient conditions for the convergence of networks composed of nodes having different 
dynamics (non-homogeneous nodes) and communicating over multiple (possibly non-homogeneous) 
media. 

3.1 The basic mathematical model and convergence analysis 

In the following, we analyze the convergent behavior of the network schematically represented 
in Figure [1] (left). In such a network, the N nodes are assumed to all share the same smooth 
dynamics and to communicate by means of the same common medium, characterized by some 
smooth dynamics: 

Xi = f{xi,z,t) i = l,...,N 

z = g{z,-^{xi,...,XN),t) 

A simplified version of the above model was recently analyzed by means of a graphical algorithm 
in (Russo & di Bernardo 2009a). In the above equation, the set of state variables of the nodes is 
Xi, while the set of the state variables of the common medium dynamics is z. Notice that the nodes 
dynamics and the medium dynamics can be of different dimensions (e.g. Xi G M", z G M"'). The 
dynamics of the nodes affect the dynamics of the common medium by means of some (coupling, or 
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input) function, \l/ : M^" — t- W^. These functions may depend only on some of the components of 
the Xi or of z (as the example in Section IGTTI illustrates) . 

The following result is a sufficient condition for convergence of all nodes trajectories of ([7]) 
towards each other. 

Theorem 4. All nodes trajectories of network ^ globally exponentially converge towards each 
other if the function f {x, v(t),t) is contracting for any v(t) G M.'^. 

Proof. The proof is based on partial contraction (Theorem [2]). Consider the following reduced order 
virtual system 

y = fiy,z,t) (8) 

Notice that now z{t) is seen as an exogenous input to the virtual system. Furthermore, substituting 
Xi to the virtual state variable y yields the dynamics of the i-th node. That is, Xi, i = 1, . . . ,N, 
are particular solutions of the virtual system. Now, if such a system is contracting, then all of its 
solutions will converge towards each other. Since the nodes state variables are particular solutions 
of ([HD, contraction of the virtual system implies that, for any i, j = 1, . . . , N: 



I cc 2 X j I y 



as t — > +00. 



The Theorem is proved by noting that by hypotheses the function f{x,v{t),t) is contracting 
for any exogenous input v{t). This in particular implies that f{y,z,t) is contracting, i.e. ([H]) is 
contracting. □ 



Remarks 

• The function (xi, . . . ,xn) is often of the form 



N 



^' (xi, . . .,xn) ■■= y^^u{xi) 



i=l 



where m : — > M'^ and all network nodes affect the medium dynamics in a similar way. 

In applications, the coupling between the nodes and the common medium is often assumed 
to be diffusive. Model ([7]) then reduces to: 



Xi = f {xi, t) + {z) - k^ {xi) i = l,...,N 

z = g{z,t) + Y,f=i [ux (xi) - (z)] 



N .... (9) 
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That is, the nodes and the common medium are coupled by means of the smooth couphng 
functions /c^ : M"^ -> M", : W ^ W and : W ^ M"^, : R'^ ^ M'^. These functions 
may depend only on some of the components of the Xi or of z (as the example in Section 16.11 
illustrates). In this case, Theorem H] implies that synchronization is attained if / {x, t) — {x) 
is contracting. Similar results are easily derived for the generalizations of the above model 
presented in what follows. 

• The result also applies to the case where the quorum signal is based not on the Xi^s themselves, 
but rather on variables deriving from the XiS through some further nonlinear dynamics. 
Consider for instance the system 

Xi = f{xi,z,t) i = l,...,N 

Ti = h{ri,Xi,z,t) i = l,...,N 

z = g{z,^ {ri,. . . ,rN) ,t) 

Theorem H] can be applied directly by describing each network node by the augmented state 
[xi, Tj), and using property (|3]) on hierarchical combinations to evaluate the contraction prop- 
erties of the augmented network dynamics; 

• Similarly, each network "node" may actually be composed of several subsystems, with each 
subsystem synchronizing with its analogs in other nodes. 



3.2 Multiple systems communicating over a common medium 

We now generalize the mathematical model analyzed in the previous Section, by allowing for s < N 
groups (or clusters) of nodes characterized by different dynamics (with possibly different dimensions) 
to communicate over the same common medium (see Figured], right panel). We will prove a sufficient 
condition for the global exponential convergence of all nodes trajectories belonging to the same group 
towards each other. This regime is called concurrent synchronization (Pham & Slotine 2007). 

The mathematical model analyzed here is 

z = g{z,-^{xi,...,XN),t) 

where: i) 7 is defined as in ([6]); ii) Xi denotes the state variables of the network nodes (nodes 
belonging to different clusters may have different dimensions, say ?T-7(i)) and z denotes the state 
variables for the common medium {z G M'^); iii) defined analogously to the previous Section, 
denotes the coupling function of the cluster 7(i) with the common medium dynamics (\I^ : x 
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Theorem 5. Concurrent synchronization is achieved in network fT^) if the functions f^(i) (x, v{t),t) 
are all contracting for any v{t) G W^. 



Proof. Recall that flTUl) is composed by nodes having dynamics /i, 
the proof of Theorem HI consider the following virtual system: 

2/1 = /i hji^z.t) 

m = f2 iy2,z,t) 

Vs = fs iVs, z,t) 



fs- Now, in analogy with 



where z{t) is seen as an exogenous input to the virtual system. Let {Xi} be the set of state variables 
belonging to the i-th cluster composing the network, and denote with Xij any element of {Xi}. 
We have that {Xij, . . . ,Xsj) are particular solutions of the virtual system. Now, contraction of 
the virtual system implies that all of its particular solutions converge towards each other, which in 
turn implies that all the elements within the same cluster {Xi} converge towards each other. Thus, 
contraction of the virtual system (fTTj) implies concurrent synchronization of the real system (fTO|) . 



To prove contraction of ffTTl) . compute its Jacobian, 



J 







dfi(yi,z,t) 
dyi 

g df2(y2,z,t) 







dy2 







... 
... 








dfs{ys,z,t) 



Now, by hypotheses, we have that all the functions fi{x,v{t),t) are contracting for any exogenous 
input. This in turn implies that the virtual system ( ITTl) is contracting, since its Jacobian matrix is 
block diagonal with diagonal blocks being contracting. □ 



3.3 Systems communicating over different media 

In the previous Section, we considered networks where some (possibly heterogeneous) nodes com- 
municate over a common medium. We now consider a distributed version of such topology, where 
each of the s < N groups composing the network have a private medium. Communication between 
the groups is then obtained by coupling only their media (see Figure [2]). The objective of this 
Section, is to provide a sufficient condition ensuring (concurrent) synchronization of such network 
topology. 
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Figure 1: A schematic representation of networks analyzed in Section \37L\ (left panel) and Section 
13.21 (right panel). The nodes are denoted with circles have a different dynamics from those indicated 
with squares. The dynamics of the common media is denoted with a rectangle. 



Note that the network topology considered here presents a layer structure. In analogy with 
the terminology used for describing the topology of the Internet and World- Wide- Web (see e.g. 
(Boccaletti et al. 2006), (Newman 2003)), we term as medium (or private) level the layer consisting 
of the nodes of the network and their corresponding (private) media; we then term as autonomous 
level, the layer of the interconnections between the media. That is, the autonomous level is an 
abstraction of the network, where its nodes' dynamics consists of the network nodes and their 
private medium. This in turn implies that in order for two nodes of the autonomous level to be 
identical they have to share: i) the same dynamics and number of nodes; ii) the same medium 
dynamics (see Figure [2]). 

In what follows we will denote with Qp the set of homogeneous nodes communicating over the 
medium Zp. We will denote with Np the set of media which are linked to the medium Zp. Each 
medium communicates with its neighbors diffusively. The mathematical model is then: 

•^i fp (-^15 ^p: -^i ^ Sp „^ 

Zp = Qp {zp, {Xp) , t) + Eje^p [<^P - (^p)] Xi^Qp 

where p = 1, . . . ,s and Xp is the stack of all the vectors Xj G Qp. We assume that the dynamical 
equations for the media have all the same dimensions (e.g. Zp e R'^), while the nodes belonging to 
different groups can have different dimensions (e.g. Xi G for any i G Qp). Here, the coupling 
functions between the media, (f)p ^ M°', are assumed to be continuous and to have a diagonal 
Jacobian matrix with diagonal elements being nonnegative and bounded. All the matrices dfp/dz 
are assumed to be bounded. 



Convergence of quorum-sensing networks 



11 



Theorem 6. Concurrent synchronization is attained in network fTB) if: i) the nodes of its au- 
tonomous level sharing the same dynamics are input equivalent; ii) fp{xi,v(t),t), gp{zp,v(t),t) are 
all contracting functions for any v{t) E M.'^; Hi) ^ are all uniformly bounded matrices. 

Proof. Consider the following 2s-dimensional virtual system, analogous to the one used for proving 
Theorem |5l 

yi,p = fp{yi,p,y2,p,t) 

y2,p = 9p (Z/2,p, Vp{t) , t) + EfeeJVp (?/2,fc) - 0p (l/2,p)] 

where p = 1, . . . , s, and Vp(t) := \E' (Xp). Notice that the above system is constructed in a similar 
way as (ITT]) . In particular, solutions of f ll2p are particular solutions of the above virtual system 
(see the proof of Theorem [5]). That is, if cluster synchronization is attained for f|T3l) . then all the 
nodes sharing the same dynamics will converge towards each other. Now, Theorem [3] implies that 
cluster synchronization is attained for system f|T3|) if: i) its nodes are contracting; ii) the coupling 
functions have a nonnegative bounded diagonal Jacobian; iii) nodes sharing the same dynamics are 
input equivalent. Since the last two conditions are satisfied by hypotheses, we have only to prove 
contraction of the (virtual) network nodes. In this view, differentiation of nodes dynamics in f[T^ 
yields the Jacobian matrix 

dfp{yi,p,y2,p,t) dfp{yi^p,y2,p,t) 

dyi,p dy2,p 
g dgp{y-L,p,Vi{t),t) 

dy2,p 

The above Jacobian has the structure of a hierarchy. Thus (see Section [2]) the virtual system is 
contracting if: 

1_ dfpiyi,p,y2,p,t) d9p{y2,pMt),t) ^^^^ contracting 
oyi,p ay2,p ° 

2. ^Myi.P,y2,p,t) bounded 
oy2,p 

The above two conditions are satisfied by hypotheses. Thus, the virtual network achieves cluster 
synchronization (Theorem [3]). This proves the Theorem. □ 

Note that Theorems H] and [5] do not make any hypotheses on the medium dynamics — synchro- 
nization (or concurrent synchronization) can be attained by the network nodes independently of 
the particular dynamics of the single medium, provided that the function / (or the /j's) is contract- 
ing. By contrast. Theorem [6] shows that the media dynamics becomes a key element for achieving 
concurrent synchronization in networks where different groups communicate over different media. 
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Figure 2: A schematic representation of the network analyzed Section [3731 The connections between 
media (and hence the connections of the autonomous level) are pointed out. Notice that only two 
nodes of the autonomous level are input equivalent since: i) their media have the same dynamics 
(in red); ii) both media are shared by the same number of nodes, sharing the same dynamics (in 
yellow) . 
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4 Synchronization control 

In Section [3], we derived several criteria ensuring node synchronization for networks where multiple 
nodes exchange their state variables using (multiple) media. The above results also allow dimen- 
sionality reduction in the analysis of the system's final behavior by treating each cluster as a single 
element, similarly to (Chung et al. 2007), a point we will further illustrate in Section |5l 

The objective of this Section is to provide a sufficient condition guaranteeing some desired 
periodic behavior for the network nodes. Specifically, we will guarantee a desired period for the 
steady state oscillations. A related problem has been recently addressed in (Russo, di Bernardo & 
Sontag n.d.), where entrainment of individual biological systems to periodic inputs was analyzed. 
We will now show the following result. 

Theorem 7. Consider the following network 

Xi / (xj, 2, t) i 1, . . . , iV A\ 

z = g{z,^{xu...,XN),t)+r{t) ^ ' 

where r [t) is a T -periodic signal. All the nodes of the network synchronize onto a periodic orbit of 
period T if: i) f {xi,v{t),t) and g {z,v{t),t) are contracting functions for any v{t) G M''; ii) ^ is 
bounded. 



Proof. Consider the following virtual system: 



yi = f{yi,y2,t) 

m = g{y2,v{t),t) + r{t) 



(15) 



where v{t) := {xi, . . . ,xn)- We will prove the Theorem by showing that such a system is con- 
tracting. Indeed, in this case, the trajectories of f|T5|) will globally exponentially converge to a 
unique T-periodic solution, implying that also Xi will exhibit a T-periodic steady state behavior. 
Differentiation of the virtual system yields: 



dfiyi,y2,t) 
dyi 




df{yi,y2,t) 

dy2 
dgiy2,v{t),t) 

dy2 



The above Jacobian has the structure of a hierarchy. Thus (see Section [2]) the virtual system is 
contracting if: 

df{y\,y2,t) dg{y2,v{t),t) both Contracting 
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2. M|L^ is bounded 

dy2 

The first condition is satisfied since, by hypotheses, the functions f {x,v(t),t) and g {z,v{t),t) are 
contracting for any v G M''. The second condition is also satisfied since we assumed df/dz to be 
bounded. The Theorem is then proved. □ 

The results can be extended to the more general case of networks of non homogeneous nodes 
communicating over non homogeneous media. 

Theorem 8. Consider the following network 



•^i fp (-^i; ^p) ^) -^i ^ Gp 

Zp = gp {zp, ^ (Xp) , t) + EfceJV„ [(l> (zk) - 4> (zp)] + r (t) Xj G Qp 



(16) 



where Xp is the stack of all the Xi G Qp and r (t) is a T -periodic signal. Concurrent synchronization 
is attained, with a steady state periodic behavior of period T if: 

1. the nodes of the autonomous level sharing the same dynamics are input equivalent; 

2. the coupling functions have bounded diagonal Jacobian with nonnegative diagonal elements; 

3. fp (xj, v(t),t) and gp {zp, v(t),t) are contracting functions for any v{t) G M'^; 
4- ^ are all uniformly bounded matrices. 

Proof. The proof is formally the same as that of Theorem [6] and Theorem [TJ and it is omitted here 
for the sake of brevity. □ 

A simple example 

Consider a simple biochemical reaction, consisting of a set of > 1 enzymes sharing the same 
substrate. We denote with Xi, . . . , X^ the concentration of the reaction products. We also assume 
that the dynamics of S is affected by some T-periodic input, r{t) (the behavior of networks where 
the medium dynamics is affected by an exogenous input will be analyzed in Section H]). We assume 
that the total concentration of Xi, i.e. Xi^T, is much less than the initial substrate concentration, 
Sq. In these hypotheses, a suitable mathematical model for the system is given by (see e.g. (Szallasi 
et al. 2006)): 

X, = -aX, + ^ z = l,...,N 

S = -EL^ + rit) ^''^ 
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with Ki and K2 be positive parameters. Thus, a suitable virtual system for the network is 



y2 z^i=i X2+J/2 



r(t) 



Differentiation of the above system yields the Jacobian matrix 

K2 







-N-, ^2 



(^^■2+^1)2 



(19) 



It is straightforward to check that the above matrix represents a contracting hierarchy. Thus, all 
the trajectories of the virtual system globally exponentially converge towards a unique T-periodic 
solution. This, in turn, implies that Xi, i = 1, . . . , N, globally exponentially converge towards each 
other and towards the same periodic solution. 

Figure [3] illustrates the behavior for = 3. Notice that, as expected from the above theoretical 
analysis, Xi, X2 and X3 synchronize onto a periodic orbit of the same period as r{t). 




10 20 30 40 50 60 70 80 90 100 




10 20 30 40 50 60 70 80 90 100 

time (minutes) 



Figure 3: Simulation of ( fTTl) . with N = 3 and r{t) = 1.1 + sin(0.1 * t). System parameters are set 
as follows: a = 1, K2 = I, Ki = 2. 
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5 Emergent properties as N increases 

In this Section, we analyze how the convergence properties of a given quorum sensing network vary 
as the number N of nodes increases. We show that for typical quorum sensing networks, as 
becomes sufficiently large, synchronization always occurs. One particular modeling context where 
these results have important implications is that of cell proliferation in biological systems. 



5.1 A lower bound on ensuring synchronization 

It is well known (Wang & Slotine 2005) that for all-to-all diffusively coupled networks of the form 



N 

Xi = f{xi,t) + ^k{xj - Xi) (20) 

i=l 



the minimum coupling gain k required for synchronization is inversely proportional to the number 
of nodes composing the network. That is, 

^min ^ 

We now show that a similar bound holds for nodes coupled by means of quorum sensing of the form 

Xi = f {xi,t) + kN{z - Xi) i = l,...,N 
z = g{z,'${xi,...,XN),t) 

To simplify notations, the above model assumes that z and all Xj have the same dimensions. Also 
note that in fl?Il) the dependence of the coupling gain on the number of nodes, A^, is given explicitly. 

Theorem 9. Assume that the Jacobian (|^) is upper-bounded by a for some matrix measure n, 
i. e., 

3a G M, Vx, Vt > 0, /i (^j < a 



Then, network [21]) synchronizes if 

a 

That is, /Cmin (xl/N . 
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Proof. Consider the virtual system 

y = f{y,t) + kN{z-y) (22) 

Synchronization is attained if the virtual system is contracting. Now, computing the matrix measure 
of the Jacobian of ( |22l) yields 



Vx, Vt > 0, fi - kN?j < fi (^^^ + kNfi (-/) < a-kN 
Thus, the virtual system is contracting if k > j^. □ 



5.2 Dependence on initial conditions 

We now consider the basic quorum sensing model ([7]). We derive simple conditions for the final 
behavior of the network to become independent of initial conditions (in the nodes and the medium) 
as becomes large. 

Theorem 10. Assume that for ^ the following conditions hold: 

• /i (1^) -oo as N ^ +00 

• g{z,V2(t),t) is contracting (for any V2(t) in W^) 

are hounded for any x, z, V2 (where \\-\\ is the operator norm) 

Then, there exists some N* such that for any N > N* all trajectories of ^ globally exponentially 
converge towards a unique synchronized solution, independent of initial conditions. 



dg 
dv2 



Proof. We know that contraction of / {x, Vi{t),t) for any fi(t) (which the first condition implies for 
N large enough) ensures network synchronization. That is, there exists a unique trajectory, Xs{t), 
such that, as t — !■ +oo, 

\xi — Xs\ ^ 0, Wi 

Therefore, the final behavior is described by the following lower-dimensional system: 

z = g{z,'^{xs),t) ^ ' 
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If in turn this reduced-order system (12311 is contracting, then its trajectories globally exponentially 
converge towards a unique solution, say xl{t), regardless of initial conditions. This will prove the 
Theorem (similar strategies are extensively discussed in (Chung et al. 2007)). 

To show that (l23l) is indeed contracting, compute its Jacobian matrix, 

dxs dz 
dg dg 
dxs dz 

Lemma [1] in the Appendix shows that the above matrix is contracting if there exists some strictly 
positive constants 6i, 62 such that 







dg 


\dxsj 




dxs 



and /i 



dg 
dz 



+ 



01 



di 



dz 



(24) 



are both uniformly negative definite. 



Now, /i ( £ 



and /i (If ) are both uniformly negative by hypotheses. Furthermore, fi ' — 



dXa 



tends to —00 as N increases: since 



1^1 

I dz I 



and 



dg 



dx. 



N* such that for any N > N* the two conditions in f l24j) are satisfied. 



are bounded, this implies that there exists some 

□ 



Also, assume that actually the dynamics / and g do not depend explicitly on time. Then, under 
the conditions of the above Theorem, the reduced system is both contracting and autonomous, 
and so it tends towards a unique equilibrium point (Lohmiller & Slotine 1998). Thus, the original 
system converges to a unique equilibrium, where all equal. 

In addition, note that when the synchronization rate and the contraction rate of the reduced 
system both increase with A^, this also increases robustness (Pham & Slotine 2007) to variability 
and disturbances. 



5.3 How synchronization protects from noise 

In this section, we discuss briefly how the synchronization mechanism provided by dynamical quo- 
rum sensing protects from noise and variability in a fashion similar to the static mechanism studied 
in (Tabareau et al. 2010). We show that the results of (Tabareau et al. 2010), to which the reader is 
referred for details about stochastic tools, extend straightforwardly to the case where the dynamics 
of the quorum variables cannot be neglected or indeed may play a central part, as studied in this 
paper. 
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Assume that the dynamics of each network element Xi in (l2Tll is subject to noise, and consider, 
similarly to (Tabareau et al. 2010), the corresponding system of individual elements in Ito form 

dx^ = {f{xi,t) + kN{z-Xi))dt + adWi i = l...N (25) 

where the all-to-all coupling in (Tabareau et al. 2010) has been replaced by a more general quorum 
sensing mechanism. The subsystems are driven by independent noise processes, and for simplicity 
the noise intensity a in the equations above is assumed to be constant. We make no assumptions 
about noise acting directly on the dynamics of the environment/quorum vector z. 

Proceeding exactly as in (Tabareau et al. 2010) yields similar results on the effect of noise. In 
particular, let x* be the center of mass of the Xi, that is 




Notice that when all the nodes are synchronized onto some common solution, say Xs{t), then, 
by definition, x* = Xs{t). 

Adding up the dynamics in fl25l) gives 

^ > ' adWi (26) 



dx' = ^ f{xi, t) j dt + kN{z - x')dt + j^Yl' 



Let 

e = /(x-,t)- 



Note that e = when all the nodes are synchronized. 

By analogy with f l25|) . equation fl26|) can then be written 

dx* = (/(x*, t) + kN{z - X*) + e) dt + (27) 

i 

Using the Taylor formula with integral remainder exactly as in (Tabareau et al. 2010) yields a bound 
on the distortion term e, as a function of the nonlinearity, the coupling gain k, and the number of 
cells N, 



E(||6||) < A^ax ( ^ ) pikN) 
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where Amax(g;^) is a uniform upper bound on the spectral radius of the Hessian and p{kN) — )■ 
as kN — )■ +00. In particular, in ( 1271) . both the distortion term e and the average noise term 
jq Si o^dWi tend to zero as — )■ +00. 

Note that an additional source of noise may be provided by the environment on the quorum 
variables themselves. We made no assumptions above about such noise, which acts directly on the 
dynamics of the environment /quorum vector z. How it specifically affects the common quantity z 
in f l25|) could be further studied. 

Similar results hold for the effects of bounded disturbances and dynamic variations. 



6 Examples 

6.1 Controlling synchronization of genetic relaxation oscillators 

We now consider the problem of synchronizing a population of genetic oscillators. Specifically, we 
consider the genetic circuit analyzed in (Kuznetsov et al. 2004) (a variant of (Kobayashi et al. 2004)), 
and schematically represented in Figure |H Such a circuit is composed of two engineered gene 
networks that have been experimentally implemented in E. coli; namely: the toggle switch (Gardner 
et al. 2000) and an intercell communication system (You et al. 2004). The toggle switch is composed 
of two transcription factors: the lac repressor, encoded by gene lad, and the temperature-sensitive 
variant of the Ac/ repressor, encoded by the gene cI857. The expressions of cI8547 and lad are 
controlled by the promoters Ptrc and Pl* respectively (for further details see (Kuznetsov et al. 2004)). 
The intercell communication system makes use of components of the quorum-sensing system from 
Vibro fischeri (see e.g. (Ng & Bassler 2009) and references therein). Such a mechanism allows 
cells to sense population density through the transcription factor LuxR, which is an activator of 
the genes expressed by the Piux promoter, when a small molecule AI binds to it. This small 
molecule, synthesized by the protein Luxl, is termed as autoinducer and it can diffuse across the 
cell membrane. 

In (Kuznetsov et al. 2004), the following dimensionless simplified model is analyzed (see Figure 

ED: 

= rr^ ^ 1+^ ~ '''' ^ 

Vi = - d2Vi (28b) 

1 + u 
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Extracellular 
autoinducer (Ale) 



Figure 4: A schematic representation of the genetic circuit: detailed circuit. 



N 



Wf. 



d^Wi I + 2d {we - Wi) 



dpWp 



Wi - We) 



(28c) 



(28d) 



i=l 



where Vi and Wi denotes the (dimensionless) concentrations of the lac repressor, A repressor 
and LuxR-AI activator respectively. The state variable We denotes instead the (dimensionless) 
concentration of the extracellular autoinducer. 

In (Kuznetsov et al. 2004), a bifurcation analysis is performed for the above model, showing 
that synchronization can be attained for some range of the biochemical parameters of the circuit. 
However, as the objective of that paper was to analyze the onset of synchronization, the problem 
of guaranteeing a desired oscillatory behavior was not addressed. In what follows, using the results 
derived in the previous sections, we address the open problem of guaranteeing a desired period for 
the steady state oscillatory behavior of network 



The control mechanism that we use here is an exogenous signal acting on the extracellular 
autoinducer concentration, see also (Russo, di Bernardo & Sontag n.d.). That is, the idea is to 
modify (128dp as follows 

Dp ^ 

" ^ (ti'i - We) - deWe + T (t) (29) 
=1 



We 



N 
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Figure 5: Simplified circuit using for deriving tlie matliematical model ( !28|) . Both the promoters 
and transcription factors are renamed. 



where r (t) is some T-periodic signal. The set up that we have in mind here is illustrated in Figure 
[6l where multiple copies of the genetic circuit of interest share the same surrounding solution, on 
which r (t) acts. From the technological viewpoint, r (t) can be implemented by controlling the 
temperature of the surrounding solution, and/or using e.g. the recently developed microfluidics 
technology (see e.g. (Beebe et al. 2002) and references therein). 

In what follows, we will use Theorem H] to find a set of biochemical parameters that ensure 
synchronization of fl28ap - (]28d|) . This, using the results of Section HI immediately implies that the 
forced network (I28ap - (]28cp . ( j29l) globally exponentially converges towards a T-periodic steady state 
behavior. 



System (128|) has the same structure as Q), with Xi = [ui,Vi,Wiy 



and: 



+ T-TT- — d,iUi 



02 

i+u7 



- d2Vi 



-d^Wp, 



kz {z) kx {xi) 








2d {We 



We know from Theorem H] that all nodes trajectories converge towards each other if: 



1. f {xi,t) — kx (xi) is contracting; 

2. g {z, t) — Nuz {z) is contracting. 
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Figure 6: Network control setup. 
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That is, contraction is ensured if there exist some matrix measures, fi^: and /x,,,*, such that 

fi^{{xi,t) - kx{xi)) and fi*^ {g {z,t) - Nu^ (z)) 

are uniformly negative definite. We use the above two conditions in order to obtain a set of 
biochemical parameters ensuring node convergence. A possible choice for the above matrix measures 
is /i* = /i** = /ii (see (Russo, di Bernardo & Slotine n.d., Russo, di Bernardo & Sontag n.d.)). 
Clearly, other choices for the matrix measures /i* and yU** can be made, leading to different algebraic 
conditions, and thus to (eventually) a different choice of biochemical parameters. 

We assume that /3 = = 7 = 2, and show how to find a set of biochemical parameters satisfying 
the above two conditions. 

Condition 1. Differentiation of — yields the Jacobian matrix (where the subscripts have 
been omitted) 



Ji :-- 



Now, by definition of /ii, we have: 



J —2a\v 2a3W 



(30) 



r 2a2U 2ea4U 2aiv 2a3W 

Hi (Ji) = max < — di H ^ H tj, —02 H 9, —£"3 — 2d-\ 

' I (1 + ^2)2 (1 + ^2)2' (l + t;2)2' (l + w;2)2 

Thus, Ji is contracting if /xi (Jj) is uniformly negative definite. That is, 

J I 2a2U I 2ea4U 
—O'l -r 77- — 5T7 "T 7;- — 577 

-^2 + (31) 

-edg -2d^ 



(IW)" 



are all uniformly negative. Notice now that the maximum of the function a (v) = ^,^^^^^2 is d = -^^f^. 
Thus, the set of inequalities ( 13T|) is fulfilled if: 



6a2\/3 I deaiVS 



-d2 + ^fi (32) 



16 ' 16 
16 

-eds - 2d + 



are all uniformly negative. 

Condition 2 In this case it is easy to check that the matrix Je := || — is contracting for 
any choice of the (positive) biochemical parameters Df.-, d^- 
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Thus, we can conclude that any choice of biochemical parameters fulfilling (1321) ensures synchro- 
nization of the network onto a periodic orbit of period T. In (Kuznetsov et al. 2004), it was shown 
that a set of parameters for which synchronization is attained is: ai = 3, 0^2 = 4.5, = 1, = 4, 
£ = 0.01, d = 2, di = d2 = d^ = 1. We now use the guidelines provided by fl32l) to make a minimal 
change of the parameters values ensuring network synchronization with steady state oscillations of 
period T. Specifically, such conditions can be satisfied by setting di = 6, ^2 = 2. Figure [7] shows 
the behavior of the network for such a choice of the parameters. 




100 200 300 400 , 500 600 700 800 900 1000 



time (minutes) 

A A ' A A' A 'A A' A I\ 



A A A 



100 200 300 400 , 500 600 700 800 900 1000 
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Figure 7: Behavior of ( l28al) -(l28cl). (l29l) . when forced by r (t) = 1 + sin (O.lt). Notice that the nodes 
have initial different conditions, and that they all converge onto a common asymptotic having the 
same period as r (t). 



6.1.1 Biological oscillators communicating over different media 

In the above Section, we assumed that all the genetic circuits shared the same surrounding solution. 
We now analyze the case where two different clusters of genetic circuits are surrounded by two 
different media. The communication between clusters is then left to some (eventually artificial) 
communication strategy between the two media (see Figure [8]). 

Notice that only one of the two media is forced by the exogenous T-periodic signal r (t) (thus 
the dynamics of the two clusters are not the same), while the two media communicate with each 
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Plate2 






^^^^^^ 




We 







Figure 8: Two clusters of genetic circuits communicating over two different media 

other in a diffusive way. The mathematical model that we analyze here is then: 

diUii 



Vil = 

Wil 

Wei 

Ui2 -- 

Vi2 = 

Wi2 

We2 



ax 



0-2 



1+u 

d: 



d2Vil 

- dsWii^ + 2d {wei - Wii) 

Wel) - deWel + r (t) + 

l+w'. 



[We2) 



[.Welj 



TT \ "n T — ^1 UiiO 



(33) 



Ol2 



1 + U 



d2Vi2 

1+^ - dzw,a ] + 2d {we2 - ^^2) 



it Ei=l (^i2 - We2) - deWe2 + (^el) " (^62) 



where xn = [ua, Vii,Wii] and Xi2 = [ui2, f i2, ^12] denote the set of state variables of the i-th oscil- 
lator of the first and second cluster respectively. Analogously, Wei and We2 denote the extracellular 
autoinducer concentration surrounding the first and second cluster of genetic circuits. In the above 
model we assume that the biochemical parameters of the two genetic circuits and media are the 
same. 

To ensure concurrent synchronization, we tune the biochemical parameters of the two clusters 
of oscillators and design the coupling function between the media (0(-)) by using the guidelines 
provided by Theorem [61 Furthermore, using Theorem [8] we can conclude that the steady state 
behavior of the two clusters is T-periodic. 



It is straightforward to check that the hypotheses of Theorem are all satisfied if: 
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the biochemical parameters of the two clusters fulfill the conditions in (13211 : 
the coupling function (■) is increasing. 



In fact, the topology of the autonomous level of the network is input equivalent by construction. 
Figure [9] shows the behavior of fl33|) when the biochemical parameters of the oscillators are tuned 
as in the previous Section, and (p (x) = Kx, with K = 0.1. 



•^60 

fH 

-+-^20 
fH 

cd 
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Figure 9: Behavior of ( 133|) when r (t) = 1 + sin (0. It). Notice that network nodes have different 
initial conditions, concurrent synchronization is attained for the network. Both the clusters exhibit 
a steady state behavior having the same period as r (t) 



6.1.2 Co-existence of multiple node dynamics 

We analyze the case where the two clusters in the previous Section are now both connected to a 
third cluster composed of Van der Pol oscillators coupled by means of a quorum-sensing mechanism. 
The three clusters have three different media, and communication between them occurs by means 
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of some coupling function. The mathematical model considered here is then: 



Uj\ — 3 h T, — 55 fli Uj-i 

I^il = ^ (ifuj; ~ dsWilj + 2d {Wel - Wil) 

^el = YliLl i^il - ^el) - deWel + (u^es) - (^el) 
7/ — "1 _L "awg, , 

- + t:;:^ - diUi2 

= if^ - d2Vi2 

Wi2 = £ (if^ - d3Wi2^ + 2d {We2 - Wi2) 

^e2 = ^f=l (^»2 - U)e2) " 4^62 + (u^es) " (We2) 
2/H = y2i 

y2i = -a {yli - 13) y2i - u^Vu + K {we^ - yu) 



(34) 



with [yu, y2i]^ denoting the state variables of the i-th Van der Pol oscillator, and with N^dp indicating 
the number of Van der Pol oscillators in the network. In the above model the Van der Pol oscillators 
are coupled by means of the medium Wes G M. The three media, i.e. Wei, We2, Wes, communicate 
by means of the coupling function We assume that the function g governing the intrinsic 

dynamics of the medium We3 is smooth with bounded derivative. The parameters for the Van der 
Pol oscillator are set as follows: a = P = u = 1. Notice that now no external inputs is applied on 
the network. 

Recall that Theorem O ensures synchronization under the following conditions: 



1. contraction of each cluster composing the network; 

2. topology of the autonomous level of the network connected and input equivalent. 

Notice that the second condition is satisfied for the network of our interest. Furthermore, 
contraction of the two clusters composed of genetic oscillators is ensured if the their biochemical 
parameters satisfy the inequalities in (132|) . 

To guarantee the convergent behavior of the cluster composed of Van der Pol oscillators, we have 
to check that there exist two matrix measures, fi^, and fi^.^,, showing contraction of the following two 
matrices: 

J = \ ° ^ 

^ -« ivi - l3)-uj^ -2ay2iyii - K 



(35a) 
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J2 



dg 

dw, 



-K 



(35b) 



e3 



Now, in (Wang & Slotine 2005), by using the Euchdean matrix measure, i.e. /i2i, it is shown 
that the matrix fl35al) is contracting ii K > a. On the other hand, to ensure contraction of J2, we 



have to choose K > G, where G is the maximum of 



da 



Thus, contraction of the cluster composed 



of Van der Pol oscillators is guaranteed if the coupling gain, K, is chosen such that: 

K > max I a, (5} 



In Figure [TOl we set g (x) = sin (x), K = 2.5, Ny^p = 2 and (x) = Kx, with K = 3. Such a 
Figure shows that concurrent synchronization of ( !34ll is attained, in agreement with the theoretical 
analysis. 
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Figure 10: Behavior of The two clusters of genetic oscillators (not directly connected) con- 

verge onto the same common evolution. Synchronization for the cluster composed of Van der Pol 
oscillators is also attained. 



6.2 Analysis of a general Quorum- Sensing pathway 

In the previous Section, we showed that our results (with appropriate choice of matrix measure) 
can be used to derive easily verifiable conditions on the biochemical parameters of the genetic 
oscillator ensuring contraction, and hence synchronization (onto a periodic orbit of desired period) 
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and concurrent synchronization. We now show that our methodology can be apphed to analyze a 
wide class of biochemical systems involved in cell-to-cell communication. 

We focus on the analysis of the pathway of the quorum sensing mechanism that uses as autoin- 
ducers, molecules from the AHL (acyl homoserine lactone) family. The quorum sensing pathway 
implemented by AHL (see Figure [TTi) is one of the most common for bacteria and drives many 
transcriptional systems regulating their basic activities. 




Figure 11: The quorum sensing pathway implemented by AHL 



We now briefly describe the pathway of our interest (see (Muller et al. 2006) for further details). 
The enzyme Luxl produces AHL at (approximately) a constant rate. AHL in turn diffuses into 
and out of the cell and forms (in the cytoplasm) a complex with the receptor LuxR. Such complex 
polymerizes and then acts as a transcription factor, by binding the DNA. This causes the increase 
of the production of Luxl, generating a positive feedback loop. 

The pathway can be described by a set of ordinary differential equations (using the law of mass 
action, see (J.Dockery & Keener 2004), (Muller et al. 2006)). Specifically, denoting with Xe the 
mass of AHL outside of the cell and with Xc the mass of AHL within the cell, we have the following 
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mathematical model: 



•^C ~^ A_^n ^C-^C d\Xf^ d2Xg 

■^thresh''' c 

Xg diX^ d2X£ 'Je-^e 



(36) 



The physical meaning of the parameters in (l36ll is given in Table [2l 

Table 2: Biochemical parameters for system (136!) 



Parameter 


physical meaning 


a 


Low production rate of AHL 


f3 


Increase of production rate of AHL 


7c 


Degradation rate of AHL in the cytosol 


7e 


Degradation rate of AHL outside the cell 


di 


Diffusion rate of the extracellular AHL 


d2 


Diffusion of the intracellular AHL 


•^thresh 


Threshold of AHL between low and increased activity 


n 


Degree of polymerization 



Now, contraction of the above system is guaranteed if 
1. —7c + -. — thresh g jg uniformly negative definite; 



2. —d2 — 7e is uniformly negative definite. 

Recall that Xc and Xe are both scalars. Now, the second condition is satisfied since system parameters 
are all positive. That is, to prove contraction we have only to guarantee that 



-7c + 



y^thresh -^c) 



is uniformly negative. Since 



-7c + 



•^thresh "^c 



thresh-^ c , 3/3 \/3 
2 ^ -7c + 



^■^thresh 

contraction is ensured if the biochemical parameters /3, g and Xthresh fulfill the following condition 

P ^ 87c 



•^thresh 3^3 
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7 Concluding remarks 

In this paper, we presented a systematic methodology to derive conditions for the global exponential 
convergence of biochemical models modeling quorum sensing systems. To illustrate the effectiveness 
of our results and to emphasize the use of our techniques in synthetic biology design, we analyzed 
a set of biochemical networks where the quorum sensing mechanism is involved as well as a typical 
pathway of the quorum sensing. In all such cases we showed that our results can be used to 
determine system parameters and dynamics ensuring convergence. 

A Proofs 

To prove Theorem [3] we need the following Lemma, which is a generalization of a result proven in 
(Russo, di Bernardo & Sontag n.d.): 

Lemma 1. Consider the block- partition for a square matrix J: 

\ A{x) B{x,y) ' 
[C{x,y) D{y) 

where A and D are square matrices of dimensions ua x and ud x rij:, respectively. Assume that 
A and B are contracting with respect to fi^ and fi^ (induced by the vector norm |«|^ and 
Then, J is contracting if there exists two positive real numbers 9i, $2 such that 

/.^(A) + |||C(x,y)|U^^<-ci 

/iM/^) + |||5(x,i/)||^_^<-4 

where \\*\\ad ^'^^ II*IIz)A ^'^^ operator norms induced by |«|^ and |«|^ on the linear operators 
C and B. Furthermore, the contraction rate is = max {c^,Cg}. 

Proof. Let z := {x, y)'^. We will show that, with the above hypotheses, J is contracting with respect 
to the matrix measure induced by the following vector norm: 

\z\ := 6*1 \x\^ + 02 \y\j;) 

with 6*1, > 0. In this norm, we have 

|(/ + hJ)z\ =9i\{I + hA)x + hBy\^ + ^2 |(/ + hD)y + hCx\^ 
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Thus, 

|(/ + hJ)z\ <9i\{I + hA)x\^ + h9i \By\^^^ + ^2 |(/ + hD)y\^ + hO^ \Cx\^^j, 

Pick now h > and a unit vector z (depending on h) such that ||(/ + hJ)z\\ = |(/ + hJ)z\. We 
have, dropping the subscripts for the norms: 



+ hJ\\ - 1) < i (\\I + hA\\ -l + ^h \\C\\] \x\ ^1 + r ( 11^ + ^^11 ■ n 
n n \ til J h \ ti2 

Since 1 = |z| = tii \x\, + 6^2 |?/| o, we finally have 



-h\\B\\^ \y\O2 



+ hJ\\ - 1) < max (\\I + hA\\ - 1 + ^/i , \ (\\I + hD\\ -l+^-^h \\B\ 

n J h \ ti2 

Taking now the limit for h 0~^: 



/i (J) < max <; fiiA)A + ^ ||C|| , + ^ \\B\ 

til ti2 



thus proving the result. 



□ 



Following the same arguments. Lemma [T] can be straightforwardly extended to the case of a real 
matrix J partitioned as 

Jii J12 . . . JlN 
J= 

JnI Jn2 ■ ■ ■ JnN 

where the diagonal blocks of J are all square matrices. Then J is contracting if 

KJu) + If ll'^12|| + • • • + I7 II^IA'II < -4l 



fJ^i-JwN) + \\Jni\\ + • • • + \\-Jin\\ < -C%j^ 

(where subscripts for matrix measures and norms have been neglected). 



(37) 



Proof of Theorem [3] 



The assumption of input equivalence for the nodes implies the existence of a linear invariant subspace 
associated to the concurrent synchronization steady state regime. We will prove convergence towards 
such a subspace, by proving that the network dynamics is contracting. Let fifhe the matrix measure 



Convergence of quorum-sensing networks 



34 



where the nodes dynamics is contracting and define: X := {xi , . . . , xj^) , F(X) as the stack of all 
intrinsic nodes dynamics, H{X) the stack of nodes coupling functions. We want to prove that there 
exist a matrix measure, /i, (which is in general different from fif) where the whole network dynamics 
is contracting. Denote with L := {lij} the Laplacian matrix (Godsil & Royle 2001) of the network 
and define the matrix L{X), whose ij-th block, Lij{X), is defined as follows: 



4 W ■■=1 



dh. 



7{i) 



dx-i 



(Notice that if all the nodes are identical and have the same dynamics and the same coupling 
functions, then L can be written in terms of the Kronecker product, ®, as (L ® with n 

denoting the dimension of the nodes and In the n x n identity matrix.) 



The Jacobian of (El) is then: 



The system is contracting if 



is uniformly negative definite. Now: 

fdF 



J :-- 



dF 

dX ^ ' 



(38) 



OF 



dF 



Notice that, by hypotheses, the matrix —L{X) has negative diagonal blocks and zero column sum. 
Thus, using ( 1371) with 6*^ = for alH, j = 1, . . . , iV, i 7^ j yields 



^^(-L{X) 







Thus: 



/i 



dF 
dX 



- L{X) < 



dF 
dX 



Since the matrix |^ is block diagonal, i.e. all of its off-diagonal elements are zero, (1371) yields: 



max < /i/ 

x,t,i 



dx 



The theorem is then proved by noticing that by hypothesis the right hand side of the above expres- 
sion is uniformly negative. 
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